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ABSTRACT 


We  present  a  new  method  for  attacking  the  radome  problem.  It  improves 
the  standard  raiome  analyses  in  the  two  aspects:  the  description  of  the 
incident  field,  and  the  curvature  of  the  radome  surface  as  explained  below, 
(a)  The  radome  is  normally  situated  in  the  near-field  zone  of  the  antenna, 
which  may  be  a  horn,  a  slot,  or  an  array.  Taking  in  to  consideration  the 
finite  antenna  size,  we  approximately  replace  it  by  an  array  of  discrete 
point  sources,  each  or  which  radiates  a  spherical  wave.  This  approximation 
is  different  from  that  used  in  conventional  techniques  in  which  the 
incident  field  from  the  antenna  is  approximately  represented  by  a  spectrum 
of  plane  waves,  instead  of  spherical  waves,  (b)  In  calculating 
the  wave  transmission  through  the  radome,  the  curvature  of  the  radome  is 
invariably  ignored  in  conventional  analyses.  Our  approach,  however, 
does  treat*  the  radome  as  a  curved  surface  b)  calculating  the  transmission 
of  a  spherical  wave  via  ray  techniques. 

In  this  part  of  the  report,  only  the  point  source  situated  inside 
the  radome  is  studied.  Extensive  numerical  results  show  that  the  curvature 
of  the  radome  may  significantly  modify  the  field  transmitted  through  the 
radome.  This  modification  cannot  be  accounted  for  by  conventional 
radome  analyses.  In  the  second  part  of  the  report,  we  wlU  consider  the 
superposition  of  point  sources  which  simulates  arrays  or  aperture  antennas 
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I .  INTRODUCTION 


Many  practical  antennas  are  covered. by  radomes,  whose  effects  on 
the  antenna  radiation  are  of  considerable  importance,  especially  in 
today's  high-performance  radar/communication  systems.  In  the  past 
quarter  of  a  century,  several  standard  analyses  have  been  devised  for 
analyzing  radome  effects.  None  of  them  is  exact,  and  improvements  are 
always  needed.  The  present  report  described  an. effort  in  this  direction. 

A  typical  radome  problem  may  De  stated  as  follows.  Let  an  aperture 
antenna  A,  for  instance  a  horn,  a  slot,  or  a  conformal  array,  radiate  a 
known  field  ^(r)  in  free  space  (see  Figure  1).  A  protective  shield  or 
radome  I  is  placed  around  antenna  A.  The  problem  is  to  determine  the 
radiation  field  £  for  the  composite  structure,  i.e.,  the  antenna  A 
radiating  in  the  presence  of  the  radome.  This  problem  has  received  a  great 
deal  of  attention  from  many  researchers  during  the  last  two  decades,  and  a 
so-called  "best  available"  method  ,for  attacking  this  problem  appears  to 
have  emerged.  A  brief  description  of  this  method  is  given  below. 

(a)  In  the  vicinity  of  I,  the  Incident  field  is  not  a  ray  field 
(locally  p.laine-wave) .  To  circumvent  this  difficulty,  let  2*  .be  resolved 
into  a  spectrum  of  plane  waves,  namely, 

E^(r)  -  dkx  dky  vJ(k)  e1^'1  (1.1) 

w(k)  -  (^]2  f  d*  f  dy  E1  (r )  e'1^ 

Here,  it  •  (k  ,k  ,k  )  is  the  direction  of  propagation  of  the  plane-wave 
x  y  z 

spectral  component.  The  spectral  wave  number  in  the  z-dlrectlon. 


(1.2) 


e\z) 


(a)  free  space 


(b)  radome  environment 


Figure  1.  Antenna  A  and  radome  Z, 


figure  2.  Two  choices  of  incident  directions:  A1  and  PI. 


(1.3) 


i 


,  /2  2  .  ,  2. 

kz  »  /u  vc  -  (k^  +  k  )  , 

may  be  real  (homogenous  plane-wave)  or  imaginary  (Inhomogenous  •'lane- 
wave).  The  weighting  factor  W(k)  is  the  amplitude  of  the  plane-wave 
spectral  component  propagating  in  the  direction  k. 

(b)  For  each  plane-wave  component,  a  transmission  coefficient  matrix 
TqOc)  for  a  flat  dielectric  slab  can  be  obtained  from  any  standard  text  on 
EM  theory.  The  subscript  zero  of  Tg  indicates  that  it  is  derived  from  the 
assumption  of  a  plane-wave  incident  field.  The  transmitted  field  £“(2)  at 
point  2  on  the  outer  surface  of  £  is  calculated  from  the  formula 


(2)  -  J  dkx  J"  dky  TQ(k)  W(k)  e1^*1 


(1.4) 


(c)  Once  E  (2)  is  known  for  all  points  on  the  outer  surface  of  Z, 
equivalent  surface  -  urrent  sources  (J(2),  K(2))  can  be  determined.  The 
convolution  of  the  source  with  the  Green's  function  gives  the  desired 
radiation  field  which  is  expressible  as 


£(r)  -  '  j  (G^  J  +  Gj  it)  da 


(1.5) 


outer  Z 


The  approach  described  above  is  of  course  theoretically , sound.  However, 
its  faithful  execution  is  impractical  because  of  the  extremely  laborious, 
numerical  integrations  in  (1.4)  and  (1.5).  In  the  well-quoted  analyses 
by,  Paris  [1]  and  Wu  and  Rydduck  [2],  the  numerical  integration  in  (1.4) 
is  avoided  by  approximating  the  transmitted  field  at  point  2  by 


tc(2 )  ;  f0(t0)  3(i0)  eik°'L' 


(1.6a) 


where  the  incident  direction  is  determined  by 


3 


a.  6b) 


k.Q  m  actual  ray  direction  Al,  or  the  direction  Pi  cf  the 
Poynting  vector  of  (Figure  2). 

Note  that  the  approximation  in  (1.6)  is  to  describe  £*  by  a  plane  wave. 
Since  the  radome  is  in  the  near  zone  of  the  antenna,  this  plane  wave 
approximation  for  described  in  [1],  [2]  dees  not  seem  to  be  a  good, one. 

In  the  present  report,  we  approach  the  radome  problem  from  a  different 
viewpoint.  Instead  of  decomposing  the  incid  nt  field  S*  into  a  plane  wave 
spectrum,  we  approximate  the  finite-sized  antenna  A  in  Figure  3a  by  an 
array  B  in  Figure  3b.  Each  element  in  array  B  radiates  a  spherical  wave. 
Those  spherical  wave  constituents,  transmitting  through  the  radome  I,  are 
superimposed  to  give  rise  to  the  desired  rad* ution  fieid  E  in  tne  far  zone. 
Thus,  the  key  step  in  the  present  approach  is  to  determine  the  transmission 
of  a  spherical  wave  through  a  curved  dielectric  shell. 

V.’e  shall  apply  geometrical  optics  to  solve  the  transmission  through 
the  curved  radome.  Results  are  presented  in  this  part  (Part  I).  In  Part 
II  of  this  report,  we  shall  discuss  how  wa  superimpose  point  sources  to 
approximate  a  f '  ite-sized  antenna  in  a  practical  radome  problem. 
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II.  DESCRIPTION  OF  PROBLEM 

The  geometry  of  the  radome  problem  under  consideration  is  sketched 
in  Figure  4.  A  point  source  at  produces  a  spherical  wave  which  goes 
through  a  curved  dielectric  shell  with  nonuniform  thickness.  Ray 
techniques  are  used  to  determine  the  field  at  point  on  a  given  surface 
outside  the  shell.  First,  let  us  describe  the  various  elements  involved 
in  the  problem. 

Coordinate  Systems  and  Time  Convention.  The  main  coordinate  system 
is  tho  rectangular  system  (x,y,z),  whose  origin  is  chosen  at  the  source 
point  Pq  and  the  z-coordinate  is  in  the  direction  of  the  beam  maximum 
of  the  antenna.  Other  coordinate  systems  at  points  P^,  p£  and  P^  along 
th«  ray  are  defined  later.  The  field  is  time-harmonic  with  the  time 
factor  exp(+Jwt)  which  is  suppressed  throughout. 

Source.  We  assume  that  tha  source  has  a  well-defined  "phase  center" 

AAA 

at  point  Pq,  tha  origin  of  the  coordinate  system  (x,y,z)t  and  radiates 
a  spherical  wave  denoted  by  (E1,  H*).  If  the  antenna  is  an  array  of  point 
sources,  it  is  necassary  to  consider  each  slement  in  the  array  separately 
and  superimpose  their  final  fields  at  tha  observation  points. 

Dlalactrlc  radome.  Tha  radome  is  a  dialactric  shall  with  nonuniform 
thickness  of  ralativa  dialactric  constant  •  c/cQ  or  rafraction  index 
n  •  * c^,  and  is  bounded  by  tha  inner  and  outer  surfaces  Z^  end  E^., 
respectively.  Tha  inner  surface  (near  the  source)  is  described  by 
the  equation: 

z  -  fj^x.y)  ,  for  <  x  <  bj^  and  Cj  <  y  <  .  (2.1) 

Tha  outar  surface  2^  it  given  by  the  aquation: 

6  . 


pspg| 


2  «  f?(x,y) 


for  a-  <  x  <  b-  and  c,  <  y  <  d. 


(2.2) 


Ic  is  not  necessary  to  know  the  analytical  fora  of  the  functions  f^(x,y) 

and  f2(x,y) .  In  computation,  only  a  set  of  discrete  data  points 

{xn,yn,fQ}  with  n  *  1,2,...,H  is  needed  for  the  description  of  f(f^  or  f,,). 

these  points  are  fitted  by  a  cubic  spline  which  gives  automatically  first 

2  2  2 

and  second  partial  derivatives  of  f,  i.e.,  3f/3x,  3f/3y,  3  f/3x  ,  3  f/3x3y, 
2  2 

and  3  f/3y  .  There  are  two  requirements  for  the  cubic-spline  fit: 

(i)  the  data  points  can  be  distributed  over  a  random  grid,  but  they  must 
be  dense  enough  to  describe  the  fine  details  of  E(£^  or  E^);  (ii)  Che 
domain  of  the  data  points  (s  <  x  <  b  and  c  <  y  <  d)  must  be  somewhat 
greater  than  the  area  of  E  in  which  the  incident  ray  is  expected  to 
intersect  the  radome. 

Observation  points.  Observation  point  P^  is  located  on  a  prespecified 
surface  £^,  which  can  be  either  one  of  the  following  two  types: 

(i)  Spherical  E^  with  center  at  Pq  and  an  Infinitely  large  radius. 

In  this  case,  is  in  the  far  field,  and  the  field  at 
calculated  by  the  ray  technique  Is  the  final  result. 

(11)  Planar  £j  which  Is  just  outside  the  radotM  end  normal  to  the 

s-axls.  In  this'  case,  we  have  to  integrate  the  field  on  Ej  to 
obtain  the  far  field. 

In  later  calculations ,  we  uee  mostly  c2„\'  spherical  E^  in  (i). 
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III.  GEOMETRICAL  OPTICS  FIELD 

For  a  given  incident  field  (2*, 2^)  generated  by  the  source  at  point 
Pq  (Fig.  4),  the  asymptotic  solution  of  the  field  at  point  is  deter¬ 
mined  using  geometrical  optics  [3],  [4].  The  method  of  solution  is 
described  belov. 

A.  Method  of  Solution 

Cousider  a  ray  in  direction  (0  ,d)  extending  from  the  source  point  Pq 
to  the  point  on  Z^.  The  source  region  (Region  'a  homogeneous  and 

A 

isotropic;  hence,  the  ray  is  a  straight  line  along  the  unit  vector  Vq^. 
First,  the  distance  rQ^  is  found  and  the  coordinates  of  point  P^  are 
determined.  Then  the  unit  vector  normal  to  the  surface  at  point  P^ 

‘  A  A  . 

is  found  (Figure  5).  The  plane  of  vectors  and  establishes  the 
incident  plane.  The  angle  between  these  two  vectors  is  the  Incident 
angle  a*.  Using  Snell's  law,  the  refraction  angle  is  obtained,  which 

A 

establishes  the  direction  of  the  transnltted  wave,  r^»  In  Region  II 
(dielectric).  The  ray  in  Region  II  la  a  straight  llna  along  Che  unit 

*  *4  *<  *  »  .  . 

vector  r12-  Three  coordinate  systems  (x^.y^.r^),  (uj.v^.N^),  and 
(xj,yj,ri2),  with  cowmen  origin  at  point  P^,  are  then  established. 

They  belong  to  the  Incident  ray,  the  surface  E^,  and  the  transmitted 
ray,  respectively. 

The  incident  field  (2^,5^)  la  split  Into  a  normally  polarized  field 
<2^® ,I^P) ,  (E-vector  normsj.  to  the  incident  plane  at  F^)  and  a  parallel 
polarized  field  (2^,2^*).  The  transritted  ftald  at  point  is  obtained 


aa  follows: 


(3.1) 


in  which  and  are  the  principal  curvatures  (inverse  ci  the  radii 
cf  curvature)  for  the  ray  pencil  in  Region  II.  They  are  found  from  the 
curvature  matrix  cf  the  transmitted  ray  at  point  P^»  as  shown  in 
Subsection  B.  The  matrix  itself  is  found  from  a  formula  involving  the 
curvature  matrix  of  the  incident  tay  and  that  of  the  surfeie  E^  at  point 
P^.  The  curvature  matrix  of  the  transmitted  ray  pencil  at  point  P^  is 
also  used  to  find  the  curvature  matrix  of  the  fay  at  point  P^  incident 
upon  the  surface  Eji 

Having  the  field  incident  upon  Z^  at  point  P2*  the  ray  direction  r^» 
and  its  curvature  matrix,  we  can  proceed,  in  a  manner  similar  to  the 
transmission  through  I^,  to  find  the  field  transmitted  through  E^  at 
(Figure  6).  Thus,  a  unit  vector  Nj,  normal  to  Z^  is  obtained,  and 
together  with  defines  the  incidence  plane  at  point  P2.  The  incidence 
angle  a*  (cos  a*  ■  N2  •  r12^  i3  t*'en  calculated.  Again,  Snell's  law  is 


invoked  to  find  the  refraction  angle  at  P2»  This  angle  specifies  the 

ray  direction  r0^  in  Region  III  (outside  the  radome).  Three  coordinate 

*  i  *  ^  *  *  A  ■'t^t* 

systems  ,  (U2»V2»N2'  and  ‘x2'y2*r23^  with  coomon  ori3in  at 

point  P2  are  then  introduced. 

The  field  (I^.H^)  incident  upon  E2  at  P2  is  resolved  into  parallel 
and  normally  polarized  fields,  from  which  the  transmitted  fields  are  found 
as  follows: 


S*» 

E?. 


n  ±ia 
t2  e2 


uCP  «  Y  r  x  ttn 

H2  Y0r23  E2 


tp  Hin 
C2  H2 


rcP 

'2 


Z0E2tl  X 


23 
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in  which 


c° _ L. 

2  1  -r  v 


Oi 


1  cos  a2 

Jn2  n  i 

cos 


C  2  “  1  +  v 


P2 


cos  a. 


v  -  -  n 

p2 


cos  a. 


(3.6) 


The  field  at  observation  point  is  then  found  from  the  transmitted  field 


at  such  that. 


?3  -  (DF^)  e  Jk°r23  , 


(3.7a) 


in  which 


DF23  “  (1  +  <»lIIr23)*1/2(1  +  (l2IIr23)”1/Z  »  (3.7b) 

and  q711  and  q^11  are  the  principal  curvatures  of  the  ray  pencil  in 
Region  III.  They  are  obtained  from  the  curvature  matrix  of  the 
transmitted  ray  at  point  P^.  This  matrix  is  obtained  from  a  formula 
already  mentioned  in  connection  with  transmission  through  For 

a  typical  factor  in  Eqs.  (3.4)  and  (3.7b),  the  following  square  root 
convention  is  used: 

'  . 

+| f | ,  if  f  is  real 

f  -  1//1  +  qr  -  •  (3.8) 

+j[f[,  if  f  is  imaginary 

It  should  be  mentioned  here  that  we  have  ignored  multiple  reflections 
in  the  dielectric  radome  throughout  out  analysis. 

B.  Details  of  the  Calculation 

Coordinates  of  the  first  refraction  point,  P^.  For  a  given  ray  leaving 
the  source  point  P.,  the  coordinates  of  the  point  Pt,  intersection  of  the 


rectilinear  ray  with  the  interface  are  given  by 


X1  *  rQl  s*n  9  cos 


y^  •  rQ^  sin  8  sin 


,  2.  2  2.1/2 
<si +  n  +  *i> 


(3.9) 


21  *  r01  -os  0 

in  which  8  and  <)  specify  the  ray  direction  in  the  spherical  coordinate 
systea  with  origin  at  Pq.  Since  point  is  on  the  surface  E^(z  *  f^(x,y>), 
we  can  write 

r^  cos  8  *  f^(r^  aia  9  cos  ^  *  r01  ?ia  9  sin  ^  ’  (3.10) 


For  given  8  and  this  nonlinear  equation  must  be  so  ved  for  r^.  Once 

A 

rQ^  is  known,  (x^.y^.z^)  are  found  from  (3.9)  and  the  unit  vector  r^  is 


*1*  *  Yly  +  *1* 


(3.11) 


Coordinate  systems  at  point  P^.  The  unit  vector  N^,  normal  to  the 
surface  at  point  P^,  is 


~  ^(*.7*2) 

H1  "  I  ^(*,7.2)1 


,  h(x,,,z)  -  *  -  £l<x,7)  ,  ?  -  Iji  +  +  |j  L 


°r  ij  •  <i  -  tu*  -  «lyj> 


(3.12) 


in  which 


®*1  2  2  1/2 
f.  -  ■—  f.  *  r— =  ,  A,  -  (1  +  ff  +’ff  )1,z 

1*  3x  ly  9y  1  lx  ly 

*1  V1 


(3.13) 


a 

(Notice  that  the  direction  of  is  chess*  pointing  away  from  the  sourca.) 

Vectors  rrt.  and  N,  specify  the  incident  plane  at  point  P, .  The  coordinate 
Ui  X  aw 


*  ^  A  A  L  A  t  A  AAA 

systems  (x*,y^,r01) ,  *xi,yl*rl2^  and  ^U1’V1’N1^  at  point  ?i  for  che 
incident  ray,  the  transmitted  ray  and  the  surface  E^,  respectively,  are 
chosen  such  that  all  three  have  one  common  coordinate  perpendicular  to  the 
incident  plane,  that  is,  v^  *  y^  *  y^  (Fig.  5).  Notice  that,  in  general, 
the  coordinate  systems  could  be  chosen  arbitrarily.  However,  the  choice 
made  here  offers  some  simplification  in  the  calculation  of  curvature 
matrices,  as  is  shown  later.  Thus, 


*  *i  *c  x  r01  .1  A 

V1  "  yl  "  yl  “  TT~ . *  L7[’(yl  +  2lfly)r  +  (X1  +  2lflx)y 

'N1  x  r0l'  1 


(3.14) 


+  (xlfly  *  ylflx)zl 

in  which  x^,  y^,  and  as^  are  given  in  (3.9)  and 

L1  "  E(yl  +  zlfly)2  +  (X1  +  Zlflx)2  +  (xlfly  *  ylflx) 2,1/2  ‘  (3*15) 

*  *•» 

Then  and  x£  are  specified  as 

“1  ■  *  *1  ■  +  !ly>  *  Vl  -  h£ly»*  +  +  4> 


+  V*1  '  *l'lx)ly  +  tyl£Xjr  +  *lElx  +  *1<4  *  i>!'> 


(3.16) 


*1  *  yl  X  r01  "  7^1  {[X1(Z1  ‘  ylfly)  +  flx(yl  +  Zl)ix  +  [yl(zl  ’  xlflx? 


♦  fly(xl  *  zl),y  "  ^zl(ylfly  +  xlflx)  +  <xi  +  y^l2* 


(3.17) 


in  which  and  are  given  in  (3.13)  and  (3.15),  respectively. 
Notice  that  both  u^  and  x^  are  in  the  plane  of  incidence. 
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Ia  order  to  specify  the  coordinate  system  for  the  transmitted  ray  at 
point  P^,  the  transmission  direction  r^  is  first  obtained.  To  this  end, 
the  incidence  and  refraction  angles  need  to  be  found.  Incidence  angle 


is  such  that 


cos  a,  «  H.  •  r_.  -  — - - — 

1  1  01  Voi 


(3,18) 


Snell's  Jaw  is  then  applied,  to  find  the  refraction  or  transmission  angle 


such  that 


■  t  -1  ,  i  -1,,  2  i.1/2 

sin  *  n  sin  ■  n  (1  -  cos  a.) 


(3.19) 


Now  r^»  like  r^,  is  in  the  plane  of  incidence  and  can  be  written  as 
r12  ■  (r12  ♦.  «1)u1  +  (r12  •  N^)N^  »  sin  o^  u^^  +  cos  N^  .  (3.20) 

A  A 

Notice  that  due  to  the  particular  choice  of  and  u^,  as  given  in 

i  t 

(3-14)  and  (3-16),  a ^  and  are  always  less  than  ir/2;  hence,  their  sines 
and  cosines  are  always  positive  real  numbers  (for  lossless  dielectric)  and 
no  sign  ambiguity  exists.  Once  r^  and  y£  are  properly  defined  (in  (3.14) 
and  (3.20)),  the  third  coordinate  direction  for  the  transmitted  ray  is 
found  from 


-  y^  x  r12  »  v^  k  (sin  aj  u^  +  cos  N^)  -  cos  a*  u^  -  sin  Ni  * 


Coordinates  of  second  refraction  point,  P2:  Having  the 
coordinates  of  point  (from  3-9)  and  (3-10))  and  the  ray  direction 
vector  in  the  dielectric  r,2,  the  coordinates  of  point  P2  can  be  found. 


(3.21) 


Thus 


r12  *  x12x  +  y12y  +  z12z  “  r12(a2x  +  S2y  +  Y2z)  *  r12r12 


*12  ’  Vl2  X2  ’  *1  +  °2r12 


y12  "  S2r12  y2  “  yl  +  S2r12 


2 12  "  Y2r12  Z2  *  Z1  +  Y2r12 


(3.22) 


Since  a 0^  and  (direction  cosines  of  r^)  are  known,  the' coordinates 
of  point  P 2  are  specified  in  terms  of  only*  On  the  other  hand,  the 
coordinates  of  point  P2  satisfy  the  equation  of  the  surface  j^2  * 

Hence , 


(zx  +  Y2ri2>  *  f2(xi  +  a2r12*  yl  +  02r12)  *  (3.23) 


This  nonlinear  equation  can  be  numerically  solved  for  r^2  and  the  result 
substituted  into  (3.22)  to  calculate  X2,  72  and  Z2. 


Coordinate  systems  at  point  PjJ  Similar  to  the  case  of  point  P^,  we 


*  *•  •»  *  “t  *t  * 

set  up  the  coordinate  systems  (*2»y2,rl2^’  ^U2,V2,N2^  an<i  (x2,y2,r23^  for 


the  incident  ray,  the  surface  l and  the  transmitted  ray,  respectively, 
with  common  origin  at  P2  (Figure  6).  Skipping  the  details,  the  results  are 
given  below: 


N„ 


Tz  (z  f2xx  ’  f2yy) 


(3.24) 


in  which 


4,  - 


<X  +  4  +  f2y>1/2 


it. 


:2x  3x 


and  f 


3f2! 


2y  3y 
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Coordinates  of  the  observation  point  P^.  Once  r y  and  the  coordinates 
of  point  P?  are  known,  the  coordinates  of  P^,  the  point  of  intersection  of  the 
ray  and  observation  surface  Ij(z  *  z^,  Zg  a  constant),  are  easily  obtained. 


*23  "  r23(a3x  +  S3y  +  V}  "  r23r 


23 


X23  *  a3r23  X3  *  x2  +  °3r23 


y23  “  63r23  y3  "  y2  +  B3r23 


Z23  "  y3r23  Z3  "  Z2  +  y3r23 


(3.32) 


but  or  the  plane  surface  Z^,  z^  *  *2  +  ^3^23  *  zg»  which  gives 


r23  “ 


z0  ”  z2 


(3.33) 


Since  cty  and  are  known,  substituting  r23  from  (3.33)  into  (3.32) 
yields  the  values  of  x y  and  Zy 

Curvature  Matrix  of  the  incident  field  at  P..  Since  the  source 

■  — —  1— ■  <— —  '  »  ...  -^..i  — ■  — .  ■  ..  .  n  ■  I  .  JL 

produces  a  spherical . field,  its  curvature  matrix  in  any  orthonormal 
coordinate  system  et  point  (here  x*  end  y*)  perpendicular  to  the  ray 
propagation  direction  rQ1  is  given  as  [3],  [4], 

fl/r. 


<  • ; 


koi 


0 

l/r 


01 


r01 


(3.34) 


in  which  r^2  is  the  radius  of  the  spherical  wavefront  at  point  P^,  and  U 
denotes  e  unit  matrix. 
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Curvature  matrix  of  the  surface  1^.  It  can  be  easily  shewn  that  the 


derivatives  of  vector  r„,  =  r,  (starting  at  the  origin  ?n  and  ending  at 

Ui  1  V 

a  point  such  as  P^  on  E^)  with  respect  to  two  independent  parameters 
and  y.  on  the  surface  E, ,  as  gl/en  below, 

i  1  1 


*rl  Sxl 


lx  3x, 


3y  3z,  » 

3^  y  +  *  •  *  *  Eus  • 


^  3r.  3x..  »  3y  „  Sz.,  « 

r  •  - —  «  p~  x  +  y  +  z  '•  y  +  f  z  , 

iy  *yx  ^ y1  ^  ly 


(3.35) 


with  z,  ■  f 


1  "  fl(xlwl) 


f lx  ’  3x, 


fly  *  3y. 


lie  in  the  tangent  plane  to  Che  surface  at  point  P^.  It  can  also  be 
easily  shown  that  the  derivatives  of  the  unit  vector  XT  normal  to  E^  at 
P^,  with  respect  to  x^  and  y^,  also  lie  in  the  sane  tangent  plane  and  are 
given  as 


Mlx  "  3^  "  ^  (-flxxX  -  flyxy) 


-f,  f.  -  f,  f, 

\  4.  lx  Ixx  ly  ly*  „ 

•  .  .2  N1 

.  A1 


*1  1 


iy  3y, 


(-f,  x  -  f,  y)  ♦ 
lxy  lyy 


-f,  f,  -  f.  f, 
lx  Ixv  Iv  Iy 


(3.36) 


in  which  is  given  in  (3.12)  as  •  — (z  *  fj^x  “  ^lyy^ 

The  letter  two  vectors  esn  be  written  lr  terms  of  the  previous  two 


vectors  defined  In  (3.35),  such  that 


Curvature  matrix  of  the  transmitted  field  at  P^.  Sow  with  the  incident 

ray  curvature  matrix  and  the  surface  curvature  matrix  known,  we  proceed  to 

tz 

find  the  transmitted  ray  curvature  matrix  Q^.  It  is  obtained  from  the 
following  equation  [4] 


k(9^)TQ^  -  k0C^)TQ&  +  h^Q  1 


hj  »  k  cos  -  k^  cos  a*  .  (3.39) 


Here  Q  \  Q*,  and  Q*  are  defied  in  (u^.v^),  (**«y*)  *nd  coordinate 

systems,  respectively.  9*  and  9^  are  matrices  which  relate  the  incident 
and  transmitted  ray  coordinates  to  the  coordinates  of  the  surface  E^  and 
for  these  coordinate  systems,  as  defined  in  Figs.  5  and  6,  they  are  given 


[X1  *  U1 


yi  *  U1 


'i 

*1  *  V1 


yl  *  Vl 


k  *  “i 


yi  •  ui 


*t 

*i  *  vi 


(3.40) 


Hotice  that  due  to  the  particular  choice  of  coordinate  systems  the  transpose, 
T  -1 

9  ,  end  the  inverse,  9  ,  of  these  matrices  are  easily  found..  Then  from 

(3.39) 

Q1  "  (9l)"ltB"leiQl0l  ♦  (co*  «1  *  °’1  co*  aJ>  Q  1HqJ)"1  .  (3.41) 

Again  it  should  be  noted  that  this  Q*  is  valid  In  (x^y^.r^)  coordinate 

system. 


•  -k.~ 


Divergence  factor  In  the  dielectric.  At  this  point  ve  can  find  DF^ 
the  divergence  factor  for  the  ray  pencil  in  going  from  point  to  point 
**  given  in  (3^4)  and  repeated  here  for  convenience  (»®«  also,  (3.8)) 

»U  -  <1  ♦  qr-12)-1/2  (1  +  Q2r12)'1/2 

in  which  and  are  the  principal  curvatures  (inverse  of  the  radii  of 

curvature)  for  the  ray  pencil  in  the  dielectric.  (Region  II)  and  are  found 
fro*  the  following  equation 

q2  -  trace  (Q*)q  +  det  (Q*)  -  0  ,  Q*  - 


or 


q  '  <qll  +  q22)q  +  (qllq22  '  q12q21>  “  0 


(3. 


ql*  q2  "  2*^11  +  q22)  -  I(qll  +  q22)  "  4(qll<l22  “  q12q2151  1  * 


Curvature  matrix  of  the  ray  pencil  lnciaent  upon  at  point  P 2 


The  curvature  matrix  of  the  ray  pencil,  in  g 
in  the  dielectric,  changes  according  to  the 


(ding  fro*  point  P^  to  P^ 
following  equation  (see  [3]), 


(Qj)’1  •  (Q^)"1  ♦  rl2 


or 


[(qJ)"1  ♦  r 


However,  this  matrix  la  valid  in  a  coordinate 


t  t 

system  parallel  to 

with  origin  at  point  Pj.  It  is  transformed  |nto  the  foliowing  matrix  in 
the  (x2,y2)  coordinate  system  (belonging  to  ihe  incident  r*y  at  point  ?,) 


(3. 


12 


tn 


-i 


Q2| 


l  * 


(B)"*1  Q* 


i  *i 
2,y2 


B 


'r7i 


(3-44) 


in  which 


B  - 


*t  *i 
*1  *  *2 

yi  •  *■> 


;c . 

i  y2 

*t  *i 
yl  *  y2 


Curvature  matrix  of  the  surface  1^,  curvature  matrix  of  the  ray  pencil 
transmitted  into  Region  III  at  point  P2,  and  divergence  factor  for  the  ray 

field  in  going  from  point  ?2  to  P^.  Following  the  steps  similar  to  those 
for  point  P1 ,  and  skipping  the  details,  we  arrive  at  the  formulas  given 
below. 


1 

”  2 

*2C2  ~  f2F2 

f2E2  ”  *2F2 

r2x’r2y 

A2 

f2G2  “  g2F2 

*2^2  ”  f2F2 

in  which 

A2  " 
E2- 


2  2 

1+  f!  -f  fi 

2x  2y 

l  +  4.  '2- 

l2xC2y 

"f2xx  f  . 

~fJzz 

,  .  llzJ CL 

&2  2 

*2 

*2  &2 

“2 

E2  -1  ^2 

«-  ‘  •  <V  V  a.  A2 

U2  \v2  r2x‘r2y 


in  which 


(3.46) 
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and  finally 


*t  *t 
x2*y2 


ej  +  (cos  a!  -  n  cos  aJ)Q.2  .  ](eb_1  (3.47) 
*i  “i  u,tv, 

*2*^2  ^  ^ 


where 


e2- 


i 

cos  a2 

0 

„t 

r  t 

COS  02 

0 

o 

1. 

• 

92“ 

o 

(3.47a) 


Divergence  factor  DF23  for  the  field  in  going  from  P2  to  P 3  is  then  given 
by  (3.7b)  and  repeated  here  for  convenience  as  (see  also,  (3.8)) 


DF23  "  (1  +  qlr23)”1/2  *  (1  +  q2r33rl/2 

in  which  and  q2  are  solutions  to  the  equation  given  in, (3.42)  with 
q^,  q32,  q,2  and  q22  being  members  of  the  curvature  matrix  of  the 


transmitted  ray: 


"t  't 
x2,y2 


11 

q12 

21 

q22 

C.  Pinal  Solution 

We  now  suHsurlze  the  final  results  obtained  so  far.  The  point 
source  at  PQ  (Figure  4)  radiates  a  spherical  wave  described  by 


“•lKLor 

E  (r,8,*)  -  ~vr-y  (P(9,*)9  +  Q(9,*H), 

-jkrr 

W  A  A 

*  YrTT^f  cos$cos9  -  Q  s±a$)x  +  (P  sin$cos8  +  Q  cos$)y  -  P  sin9z] 
0  (3.48a) 


•+i  '  -*t 

-  YQr  *  E 


(3.48b) 


where  (r,0,$)  are  spherical  coordinates  with  origin  at  Pq.  The  pattern 
functions  P(0,$)  and  Q (9 , )  in  (4.28)  are  given.  At  point  (Figure  5), 

»i 

we  decompose  the  field  into  two  components  in  the  directions  of  (x. ,yj) ,  i.e., 

-  I1(r1.01,*1)  -  (Ej-ij)ij  +  (E*-y*)y*  ,  ^  -  YQr01  x  g*.  (3.49) 

"i  *i 

where  (x^.y^)  are  defined  in  (3.14)  and  (3.17).  At  the  observation  point 
P3,  we  express  the  field  as  follows 

^  .  sj  -  v23 .  q  .  (3.50) 

*£  Afc 

where  (Xjtyj)  are  defined  in  (3.25)  and  (3.31).  The  two  components 
of  in  (3.50)  are  found  from  the  matrix  equation 


.'t 

E3*X2 


(DF12)(DF23)e 


rjk0(nrt2^23) 


^P.P /-'i  “t.  n  p/i  ~t, 

'WV*!5  nt1t|(x2*y1) 

-1  p  n,~l.  At.  n  n,‘i  At. 
n  t]Lt2^Y2*xl^  ele2'y2*yl' 


li'il 
E1  X1 


;i  “i 

Vyi 


(3.51a) 


or  more  compactly, 

-jk-(nr12  +  r2_) 

1?3  -  (DF)  e  0  12  23  T  f 


(3.51b) 
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In  (3.51),  a  is  the  refraction  index  of  the  dielectric,  kg  is  the 
free-space  wave  number  aud  tj,  t^,  t^,  t^  are  the  normal  and  parallel 
transmission  coefficients  at  point  and  P^,  respectively,  as  given 
in  (3.2)  and  (3.6),  The  two  divergence  factors  are  given  in  (3.4)  and 
(3.7b).  Their  calculations  constitute  the  major  effort  of  the  present 


solution 


IV.  FAR  FIELD 


For  Che  radome  problem  under  consideration  (Figure  4) ,  ve  are 
generally  interested  in  the  transmitted  field  in  Che  far  zone  (oucside 
Che  radome),  namely, 

-jkor  . 

^(r)  "  JtfxJ  [9Pt(9^)  +  *  Qt(e,*)]  ,  r  .  (4.1) 

0  y 

where  (r,0,$)  are  Che  spherical  coordinates  with  origin  at  (Figure  4). 
From  the  analysis  in  Section  3,  we  have  found  the  field  at  an  observation 
point  P.j  over  a  surface  Z^.  Now,  let  us  consider  how  to  obtain  the  far 
field  lT"  from  the  field  over  E^.  There  are  two  ways  for  doing  this. 

A.  Direct  Ray  Method. 

Referring  to  Figure  7,  we  choose  the  surface  E^  to  be  a  sphere 
centered  at  PQ  and  with  radius  r,  where  r  -*■'  ®.  Then  point  P^  is  already 
in  the  far  zone.  The  field  at  P^  is  calculated  from  (3.51).  Because 
of  the  fact  that  P^  is  at  an  infinitely  large  distance  from  the  soiree, 
we  can  use  some  approximations  for  distance  r23.  As  shown  in  Figure  7, 
ray  A  is  the  actual  ray  (which  follows  the  Snell's  Law)  going  from 
PQ  to  P^,  whereas  ray  A'  is  simply  a  straight  line  connecting  Pq  and 
P^.  When  appears  in  the  phase  calculation  of  a  field,  we  use 


t23  '  r2'3  ”  r03  “  r02'  “  r  "  r02' 


(4.2a) 


When  r^  appears  in  the  amplitude  calculation  of  a  field,  we  use 


r23:  r  (4.2b) 

where  r  is  the  radius  of  sphere  Z^  and  is  infinitely  large.  For  example, 
we  use  (4.2b)  in  (3.7b),  and  obtain 


x  •+■  ® 


(4.3) 
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DF23  7fql  q2  ) 


This  method  for  calculating  the  far  field  is  the  simplest.  However,  it 
fails  if  the  observation  point  is  a  caustic  point  of  the  transmitted 
ray  field.  An  example  is  shown  in  Figure  3.  The  incident  field  on  the 
radome  is  a  set  of  parallel  rays  from  a  parabolic  reflector.  After  going 
through  a  dielectric  slab  radome,  the  transmitted  rays  are  still  parallel. 
Hence,  q***  *  q^j**  3n<*  *^23  *n  ^*3)  becomes  Indeterminate*.  The 

failure  of  the  direct  ray  method  does  not  occur  often  in  practical 
applications.  When  it  does,  we  may  use  the  second  method  described  next. 
B.  Fourier  Transform  Method 

Instead  of  going  to  the  far  field  directly,  we  first  calculate 
the  field  over  a  planar  surface  E^  just  outside  the  radome  by  the 
ray  method  (Fig.  9).  It  may  be  shown  that  the  far  field  EC  In  (4.1) 
is  related  to  the  Fourier  transform  of  5^.  The  exact  relation  is 
stated  below. 


Pt(9,1^)  -  At  e 


+jk0z0c°s9 


(f  cos  $  +  f  sin  $) 
x  y 


.t,fl  . .  1  cos  0  +^kOzOc030 

)  (9,$)  -  J — s -  e 


(fy  cos  4  -  fx  sin  <fr) 


(4.4a) 

(4.4b) 


where  zn  is  the  distance  from  P.  to  plane  E,.  The  two  functions  f 
0  0  3  x 

and  f  are  Fourier  transforms  of  the  field  over  namely. 


<9,$)  -  |  | 


Ex(x,y)  e 

y 


jkQ(x  sin  9  cos  *  +  y  sin  0  sin  *)  ^  (4>5) 
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where  (E^.E^)  are  tangential  components  of  .  s  given  m  (3.50).  The 
Zq  cos  0  term  in  (4.4)  is  the  correction  made  for  the  transfer  of  the 
field  with  respect  to  the  origin  at  C  to  the  field  with  respect  to  the 
origin  at  Pq  (Figure  9) . 

The  integration  in  (4.3)  is  over  the  infinitely  large  surface  of 
It  is  obvious  that  in  actual  computation  the  integration  area  Lost  be 
limited  to  the  region  where  the  fields  are  most  significant.  To  appreciate 
this  better,  we  rewrite  it  as 


fx(u,v)  -  A*  [  [  Ex(x',y')  ej2ir(x,u47’v)dx*dy* 

XT  •  XT 


(4.6) 


x'y’ 


with 


x’  ■  x/An,  y’  »  y/A  ,  u  *  sin  0  cos  <p,  v  *  sin  0  sin  <f>  .  (4.7) 


In  effect,  the  function  f  (u,v)  is  the  plane-wave  spec  tram  of  the  field 

y 

distribution  Ex(x',y')  on  the  plane  E This  field  distribution  can  be 

y 

considered  to  be  essentially  band-limited,  i.e.,  the  energy  is 

concentrated  in  the  visible  region  of  the  spectrum.  Therefore,  f  has 

’  ,  y, 

significant  values  for  ]u|  _<  1  and  |v|  <_  1. 

The  integration  in  (4.5)  is,  best  computed' by  the  use  of  the  Fast 

Fourier  Transform  algorithm.  But,  for  this,  a  uniform  rectangular 

samoling  grid  is  xequired.  So  the  first  step  would  be  to  find  the  field 

on  the  plane  E^  at  a  set  of  uniformly  spaced  points  from  the  fields  given 

at  a  set  of  randomly  spaced  (in  the  general  case)  points  which  are  the 

result  of  ray  tracing  through  the  radome  (fields  at  points  like  P^) . 

This  requires  interpolation.  We  have  developed  a  rather  efficient  linear 

interpolation  algorithm  in  the  general  two-dimensional  case.  The  grid 

constants  (distance  between  two  consecutive  points)  in  the  x  and  y 

directions  are  determined  by  the  Nyqulst  Sampling  Theorem.  Thus, 
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Therefore,  the  natural  choice  for  the  integration  points  would  be  a  square 
unifora  saapling  grid  with  grid  conatant  of  the  order  of  X^/ 2.  However, 
in  practice,  since  the  beam  angle  without  redone  is  known  and  is  not 
expected  to  change  drastically  in  the  presence  of  the  radoae,  a  better 
estlnate  of  the  upper  liaits  of  | u f  •  -j  and  jv|  »  ~  can  be  made.  These  will 
be  generally  saaller  chan  unitv  as  given  in  (4.8)  and  (4.9),  which  will 
allow  for  larger  saapling  intervals  Ax  and  Ay. 

A  reasonable  estimate  of  the  grid  size  X,  Y,  in  the  x-  and  y-d loans Ions 
can  be  made  on  the  basis  of  the  significance  of  field  values  on  plane  Z ^ 
such  that  the  field  has  appreciable  values  for 

"  f  —  *  <  2  »  *  (**W) 

X  and  Y,  together  with  Ax,  Ay,  specify  Che  number  of  grid  points  in  the 
x  and  y  directions, 

H-jf  ,  H-jI  (4.11) 

which  aust  be  the  ssae  as  ths  nuaber  of  points  in  the  u-  and  v  directions 
In  tha  transfor*  doaain.  Since  Cha  grid  axtant  in  the  transform  dosain 
0  snd  V, srs  slrasdy  spacified,  tb.*  saapling  ictarval  in  tha  transform  dome,; a 


•  ¥ 

would  ba 

0 

Au-g  , 

Av  •  — 

N  * 

(4.12) 

which  complacaa  tha 

picture. 

■ 
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To  summarize  the  main  points  of  this  section,  starting  from  a 
knowledge  of  the  field  on  a  set  of  generally  random  points  on  the 
plane,  an  interpolation  scheme  is  used  to  obtain  the  field  values  on 
M  *  H  uniformly  spaced  points  in  the  x  and  y  directions  (see  (4.11)) 
vith  the  grid  size  X  *  Y  centered  at  the  origin  and  sampling  intervals 
Ax,  Ay  in  the  x  and  y  directions  (see  (4.8),  (4.9)),  respectively.  Once 
this  is  accomplished,  an  FFT  algorithm  is  used  to  find  the  far  field. 
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V.  SPECIAL  CASES 

The  solution  obtained  for  the  present  radome  problem  is  the  most 
general  one.  Let  us  now  concentrate  on  a  few  special  cases  which  bring 
out  interesting  physical  phenomena. 

A.  Sc  radoae. 

In  the  absence  of  a  radoae  (n  »  1),  the  rays  in  Figure  4  becomes 
straight  lines  going  from  PQ  to  P^,  and  che  field  at  P^,  according  to 
(3.51),  becomes 


±  ,n„,  -Jk0(r12+r23>  * 

E,  -  (DF)  e  Ex  . 


(5.1) 


Here  the  total  divergence  factor  is 


01 


(5.2) 


DF  -  (OF  )  (DP  )  -  r—  +  r  +  r 

ul  rl2  r23 

which  accounts  for  the  spherical  spread  of  the  incident  field. 

B .  Dielectric  Slab  Radome 

When  1^  and  are  parallel  planar  surfaces,  the  radoma  shell  becomes 
a  dielectric  slab  (Figure  10).  The  field  at  is  given  by  (3.51)  with 


DF  - 


01 


t<r01  *  n  r12  *  r23)(r01  +  n6  r12  +  r23Jl1  2 


(5.3) 


where 


cos  c. 


cos  a. 


.2  i 


-  (sec  aj)[l  -  A*  ,in‘  aJ] 


12  _,_2  i. 


(5.4) 


Details  of  the  derivation  of  (5.3)  are  given  in  Appendix  A.  When 


•  0  (norael  incidence),  (5.3)  becomes 
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DF 


"01 


L01 


(l/r)r12  +  r23 


(5.5) 


It  Is  veil-known  that  an  interpretation  of  OF  i3  the  expansion  ratio 
defined  by 


DF  - 


Area  at  P. 


Area  at  P, 


1/2 


(.5.6) 


where  "Area"  is  that  of  a  small  ray  tube  centered  around  the  ray  shown  in 
Figure  10.  As  a  numerical  example,  let  the  thickness  of  the  slab  be 
1  meter  and  both  and  P^  be  1  mater  away  from  the  slab  (Figure  10). 

In  the  absence  of  the  radome  (n  ■  1),  DF  in  (5.2)  has  the  numerical 
value 

(DF)n-l*3  *  for  all  .  (5.7) 


When  the  slab  radome  is  present,  we  plot  DF  in  (5.3)  as  a  function  of 
n  for  the  two  values  of  a*  -  0*  and  45*  (Figure  11).  We  note  that  for 
n  >  1,  DF  is  always  higher  than  the  free-space  value.  In  the  limit 
n  ■*  •,  DF.  approaches  the  asymptotic  value  of  0.5  for  all  Incident 
angles,  a*. 

C.  Spherical  Shell  Radome 

Let  the  two  surfaces  and  Ej  be  spheres  of  radii  R^  and  Rj, 
respectively  (Figure  12).  The  thickness  of  the  radome  along  the  z-axis 
is  d.  Ve  concentrate  on  the  far  field  along  the  z-axis,  i.a., 

'23**  •  al*°  '  •  '  <5-8> 

The  field  at  Pj  is  given  by  (3.51b)  where  the  divergence  lector  is 


reduced  to 
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(5.9) 

tfe  note  that  the  first  factor  in  (5.9)  is  identical  to  the  DF  in  (5.5) 
for  a  dielectric  slab  radome.  Thus,  the  ratio  of  the  electric  field  at 


P^  for  the  spherical  shell  and  that  for  a  dielectric  slab  is 


|  E j  for  shell  i 

,  +  .  ,v  f0lf 

r*2-*i 

n  *  .♦  ,  " 

| 2- |  for  slab 

!  +  (n  1)  ^ 

a 

-  1  + 


1  - 


*1 


01 


-1 

l  .(5.10) 


As  a  numerical  example,  consider  the  case  in  which  E^  end  Z^  are  concentric 
(R^  -  R^  “  d),  and  (R^/d)  »  2.  Ve  plot  n  as  •  function  of  (r^/d)  for 
n  -  0.5  and  n  -  3  (Figure  13).  Not';  that  n  *  1  when  rQ;J  -  R^.  (All 
three  points  A^,  A^,  Pq  in  Figure  10  become  one  point.)  Then  the 
transmitted  field  through  a  concentric  spherical  shell  and  that  through 
a  slab  oecome  the  same. 

A  most  interesting  phenomenon  occurs  when  four  parameters 
(n,  R^,R2,r0.)  satisfy  the  following  relation 


R 


(°  -  1)coA  f.  - 1 

2  RI  +  (a  -  l)rQ1  (  a 


(5.11) 


Then  DF  in  (5.9)  becomes  infinite!  It  means  that  the  rays  in  the  pencil 
near  the  axis  emerged  from  the  radome  are  parallel  so  that  they  focus  at 
the  far  field  point.  According  to  the  present  geometrical  optics  theory, 
under  the  condition  in  (3.11),  the  far  field  on  the  z-axls  is  Infinitely 
large  (a  caustic  point  of  the  geometrical  optics  field) .  The  actual  field 
is  large  but  finite,  end  its  value  can  be  predicted  only  by  using  a  more 
refined  theory  than  the  present  geometrical  optics  (such  ar  the  Fourier 
transform  method  in  Section  4.B).  This  subject  will  be  investigated 
in  Part  II  cf  this  report. 
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Figure  13.  Transmitted  field  S  transmitted  through  e  spherical  shell 
normalized  by  which  is  that  through  a  dielectric  slab. 


VI.  NUMERICAL  RESULTS  OF  RADOME 


For  a  finite-sized  antenna  inside  a  radome,  on**  approach  is  to 
represent  the  antenna  by  an  array  of  point  sources.  In  the  present  part 
(Fart  I),  we  consider  only  a  single  point  source,  while  the  array  will  be 
studied  in  Part  II.  As  explained  in  (3.48),  the  incident  field  from 
the  point  source  is  characterized  by  two  pattern  functions  F  and  Q.  For 
the  present  computations,  we  assume  that  the  point  source  is  y-polarlzed. 
Then  it  follows  that 


P(9,*)  -  Vg(0)  sin  *  (6.1a) 

Q(0,*)  -  Vh(9)  cos  *  (6.1b) 

where  and  Vfl(0)  are,  respectively,  the  E  and  H  plane  patterns,  and  they 
assume  the  form 


Vg(9)  -  (cos  0)‘ 


vH(0)  -  (cos  0) 


(6.2) 


In  particular  for  a  -  1  and  n  -  0,  the  incident  field  is  identical  to 
the  far  field  of  a  y-directed  electric  dipole.  In  the  E-plane  ($  ■  90°), 
the  Incident  field  in  (3.48)  becomes 


-JV 


l^r.e,*  -  5v’)  «  777x"7  cos  ®1 


(6.3) 


The  total  transmitted  field  through  the  radome  in  the  E-plane  is 


given  by 


-Jk0r 


^(r.d,*  -  90*)  *  Tr'/lq)'  *®  PC0)1  •  r  -  * 


(6.4) 


In  the  following,  we  present  results  of  P(0)  for  various  radomes. 


We  have  studied  a  total  of  8  radomes.  They  differ  in  the  following 
parameters,  as  listed  in  Table  I. 

(i)  Radome  shape  (sphere  or  paraboloid) 

(ii)  Relative  dielectric  constant  er  (from  2.-5  to  5) 

(iii)  Radome  thickness 

(iv)  Source  positions 
We  discuss  our  numerical  results  belov. 

All  the  spherical  radomes  (A  to  D)  have  a  radius  of  20  Xq  (Figure  14) . 

We  plot  the  magnitude  of  pattern  function  P(0)  defined  in  (6.4)  as  a  function 
of  0  in  Figures  15  to  22.  Generally,  the  effect  of  the  radome  is  predicable, 
namely,  (i)  |P(0)|  decreases  for  tnicker  layers  or/and  higher  (ii)  |P(0)| 
decreases  as  the  source  moves  laterally  away  from  the  z  axis.  The  effect  of 
the  radome  curvature  can  be  seen  from  Figure  15.  When  the  source  is  located 
at  position  3  (center  of  the  spherical  shell) ,  the  field  on  the  z-axis  is 
identical  to  that  of  a  dielectric  slab  radome,  as  can  be  predicted  from  (5.10) 
This  field  becomes  stronger  than  its  counterpart  through  a  slab  for  the  source 
position  4,  and  becomes  weaker  for  the  source  position  2. 

Two  types  of  paraboloidal  radomes  are  considered.  In  both  types,  the 
inner  surface  is  described  by 

ff-]  -  50  -  -iy  (x2  +  y2)  (6.5) 

* 0J  8Xg 

The  outer  surface  for  the  first  type  ia 

f-1-  50.5  -  -ir  (x2  +  y2)  ,  (6.6) 

lX0J  9XJ 

so  that  the  thickness  of  the  radoce  increases  toward  its  base  (Figure  23). 

The  outer  surface  for  the  second  type  is 


TABLE  I 


RADOME  PARAMETERS 


Radome 

Shape 

er  ' 

Thickness 
along  j-axis 

Patterns  in 
Figures 

A 

15,  16 

B 

sphere 

l  .  -> 

V4 

17,  18 

C 

5  0 

V2 

19,  20 

D 

V4 

21,  22 

E 

Paraboloid  I 
(Figure  23) 

V2 

25,  26 

F 

Paraboloid  II 

2*5 

V2 

27,  28 

G 

(Figure  24) 

V4 

.  31,  32 

H 

Paraboloid  I 
(Figure  23) 

5.0 

It 

V2 

33,  34 

Aq  is  free-space  wavelength 


|P(0)| 


POLAR  ANGLE  ,  8  (OEG) 


Figure  17.  E- plane  rediet Ion  pattern  through  redone  B  (cf  »  2.5,  d 


46 


I  P(0)  I 


POLAR  ANGLE,  9  (DEG) 


Figura  18.  E-plana  radiation  pat tarn  through  radoaa  8  (cf  •  2.5,  d  •  XQ/*) 


IP(0)I 


lp(0)l 


Figure  20.  E-plane  radiation  pattern  thiough  redone  C  (cf  •  5.0,  d  •  XQ/2). 


POLAR  ANGLE  ,0  (DEG) 


’igura  22.  E-plitaa  radiation  patcarn  through  radoua  0  ( 
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-  50.5  - 


8.321A: 


,  2  .  2. 

(x  +  y  ) 


(6.7) 


so  that  the  radome  thickness  is  nearly  uniform  (Figure  2A).  For  the 
incident  field  in  (6.3),  the  E-plane  patter  |P(9}[  has  been  calculated 
for  five  source  locations  incide  radomes  E  to  H. 

It  is  particularly  interesting  to  note  that,  for  the  field  on  the 
z-axis,  there  is  a  striking  difference  between  radome  E  and  radome  F, 
namely,  there  is  a  dip  in  Figures  25  and  26  for  radome  E,  whereas  there 
is  a  peak  in  Figure  27  and  28  for  radome  F.  This  is  so  despite  the  fact 
that,  near  the  tip,  radomes  E  and  F  are  quite  similar  as  may  be  seen  from 
the  expanded  graph  in  Figure  29.  The  reason  for  such  an  anaaulous  behavior 
is  explained  below.  On  the  z-axis,  the  geometrical  parameters  of  radomes 
E  and  F  are  listed  in  Table  II.  R^  and  are  the  radii  of  curvature  of 
the  inner  and  outer  radome  surfaces. 


TABLE  II 
FIELD  ON  z-AXIS 


Radome 

VX0 

roi/x 

|F(8)| 

E 

A.  500 

50 

0.647 

F 

4 

4.161 

1.036 

respectively.  Note  that,  while  R^  is  the  same  for  both  radomes,  Rj  does 
vary  slightly  (about  8%).  According  to  (5.11),  the  critical  value  of  Rj 


for  source  position  1  (r^  -  50  A^)  is 


Critical  Rj  -  3.7  A 


At  this  critical  R2,  the  divergence  factor  DF  and,  therefore,  field 


(6.8) 


|P(0)j 


POLAR  ANGLE ,  9  (OEG) 

Figure  25.  E-plane  radiation  pattern  through  radome  E  (e  -  2.5, 
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radiation  pattern  through  radorae 


radiation  pattern  through  radoae 
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in  (5.1),  as  predicated  by  the  present  geometrical  optics  theory  becomes 
infinite.  In  Figure  30,  we  plot  the  ratio  of  the  field  with  or  without 
the  radome  for  three  values  of  r^.  the  solid  curve  iu  Figure  30  corresponds 
to , the  case  discussed  in  Table  II .  For  both  radoaes  E  and  F,  their  values 
of  are  close  to  the  critical  value  “  3.7  X  and  consequently,  the 

fields  on  the  z-axis  are  quite  sensitive  to  !L.  For  radoae  E,  the  pattern 
function  |pj  has  the  value  0.647  (dip),  whereas  that  of  radoae  F  has  the 
value  1.036  (peak).  It  should  be  pointed  but  that  the  value  1.036  for  radome 
F  is  not  a  very  large  number.  We  do  not  have  a  caustic  in  the  far  field. 
Hence,  our  calculations  near  the  peaks  in  Figures  27  and  23  based  on  the 
geometrical  optics  should  be  reasonably  accurate. 

Patterns  for  radomes  G  and  H  presented  in  Figures  31  to  34  exhibit 
the  similar  peak  and  dip  phenomenon. 
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Figure  30.  Variation  of  axial  field  atrength  with  for  a  paraboloid  with 
the  aource  on  r-axi8. 
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Figure  31.  E-plene  radiation  pattern  through  redoes  G  (e^  -  2.5,  d  -  XQM) . 
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APPENDIX  A 


DIVERGENCE  FACTOR  FOR  A  DIELECTRIC  SLAB  RADCME 

Figure  10  shows  the  slab  radoae  in  Cartesian  coordinates.  Since  the 
curvature  of  the  radoae  is  zero  everywhere,  the  curvature  matrix  of  the  radome 
is  zero.  Considering  the  transmission  at  P^,  if  is  the  curvature  matrix  of 
the  incident  wavefront,  the  curvature  matrix  of  the  transmitted  wavefront 
at  Px  is  (3.41) 

Q1  “  K  Qi  ei/nJ  (8i)”1  (A.D 

where  0*  and  0*  are  given  by  (3.40). 
i  i  ^ 

The  factor  [8^  a£/n]  may  be  evaluated  using  (3.41)  and  (3.34)  to  yield 


Bence,  (A.l)  simplifies  to 


£  i  2 

where  8  *  [cos  a^/cos  a^j  .  Nov,  (q^,  q^)  the  principal  curvatures  of  the 


transmitted  wavefront  at  are  (from  3.42) 


qX  "  l/(enr01) 
q2  "  1/(nroi) 


(A.  3) 


The  divergence  factor  from  to  i,*.,  DF^  is»  therefore,  given  by 
r  \2  1  ,  /m 


B(or01)^ 

?12  " 

|_(8nr01  +  r12)(nr01  +  r12) 


(A. 4) 


Now  let  us  consider  the  transmission  at  P2 .  The  curvature  matrix  Q2  of 
the  incident  wavefront  at  P2  is  given  by  (3.44)  where  B  is  the  coordinate 
transformation  matrix  from  (x^,  y*)  to  (x2t  y2) .  However,  in  the  present 
case  B  turns  out  to  be  a  unit  matrix  and  hence 


•413  >  %  •  -"to* .  ;i  . 

Therefor.,  _  ;i  -  [(qJ)'1  *  r^l'1 

(Q^)*1  can  be  obtained  from  (A. 2)  as 


(A.  5) 


(Qj)"1 


(A.  6) 
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Hence, 


<4  xl  ^  “ 

1  XV  y2 


(8nr0i  +  r12) 


(A.  7) 


“01  +  r12 


Again,  since  the  curvature  matrix  of  is  =,jro»  the  curvature  matrix  of 
the  transmitted  wavefront  at  P2  is  given  by  (3.47) 


(A.  8) 


where  9*  and  9^  are  given  by  (3.47a).  Equation  (A.8)  leads  to  the  following 
expression  for  oHht  ~t  : 

*  I *2  *  y 2 


.Ate 

2 


(6nr0i  +  r12) 


(A.  9) 


nrQ1  +  r1? 


Now,  as  before,  we  calculate  the  principal  curvatures  of  the  transmitted 
wavefront  at  Pj.  Thus, 


1  8nr01  +  r12 


*  '“r01  +  c12> 


(A. 10) 


And,  the  divergence  factor  DF^  in  passing  from  to  P^  is  given  by 


r  .  1  (8nr01  *  r12)(nr01  *  ,r12> 

23  "  (gnr01  +  r12  +  Bnr^) (nrQ1  +  r12  +  nr23) 


(A. 11) 
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